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Self- Consistent Fluctuation Theory for Strongly Correlated Electron Systems 

Hiroaki Kusunose* 
Department of Physics, Ehime University, Matsuyama, Ehime 790-8577 

A self-consistent theory for two-particle fluctuations with renormalized irreducible vertices is proposed. 
Using the Parquet formalism, we construct the fully antisymmetric full vertex in terms of the two-particle 
fluctuations in the charge, the spin and the particle-particle channels on an equal footing to satisfy the 
Pauli principle. The fluctuations are determined self-consistently, which are reflected into the one-particle 
self-energy via the Schwinger-Dyson equation. We demonstrate the application of the present theory to the 
impurity Anderson model and the Hubbard model on a square lattice mainly for the particle-hole symmet- 
ric case. In both models the vertex renormalization in the spin channel eliminates magnetic instabilities 
of mean-field theory to ensure the Mermin- Wagner theorem. The present theory gives the same critical 
exponents of the self-consistent renormalization theory in the quantum critical region. 

KEYWORDS: two-particle self-consistency, parquet equation, self-consistent renormalization theory, impu- 
rity Anderson model, Hubbard model, Mermin-Wagner theorem, quantum critical point 



1. Introduction 

A variety of phenomena such as a metal-insulator transi- 
tion'^^ and an anisotropic superconductivity^"'"' emerges by 
strong electron-electron correlations. The fundamental mod- 
els like the Hubbard and the Anderson lattice models have 
been investigated extensively using various analytical and nu- 
merical techniques. It has revealed that microscopic detail of 
the non-interacting band structures is capable to yield various 
type of metallic states and superconducting pairings in spite 
of simplicity of the models. Insulating phases with or without 
a magnetic ordering are described by the models as well. 

The attempts to describe effects of strong correlations in 
metallic phases often begin with the Hartree-Fock (HF) and 
the random-phase approximation (RPA). Then, the low-order 
perturbation theory and the fluctuation exchange (FLEX) ap- 
proximation""'-'' that includes systematically the collective 
fluctuations improve the treatment of correlations. They are 
mainly based on the approximation for the one-particle self- 
energy, and the latter follows the Baym and Kadanoff 's con- 
serving arguments''* '^' within the one-particle level. '^' So- 
called the dynamical mean-field theory (DMFT) taking ac- 
count of local fluctuations exactly'^"'^' is also classified into 
the one-particle based theories. The theories mentioned above 
have achieved considerable success in some aspects, however, 
at the same time they have significant inconsistency in the 
two-particle level such as a violation of the Mermin-Wagner 
theorem.^"'^'' This is inherent from the one-particle based 
theories and/or the lack of (spatial) vertex corrections.^^' ^■'' 
These theories also violate the Pauli principle that is the es- 
sential ingredient behind the Fermi-liquid (FL) theory. ^''"^^' 
We also note that an appropriate description for a magnetic 
instability especially near a quantum critical point (QCP) fa- 
cilitates a better understanding of anisotropic superconductiv- 
ity mostly appeared near the QCP.^'^' 

Meanwhile, the theories based on the two-particle self- 
consistency have been developed. The successful phe- 
nomenological method is the self-consistent renormalization 
(SCR) theory by Moriya and co-workers.^- ^■^^-3°' The SCR is 
thermodynamically consistent due to the self-consistent treat- 



ment of spin fluctuations in the mode-mode coupling term, 
and it gives the non-trivial exponents near the QCP. Since 
the SCR theory relies on the phenomenological expression 
for the spin fluctuation that is valid in the FL regime, a 
range of the applicability of the SCR is a delicate issue espe- 
cially in heavy-fermion systems where a crossover to a local- 
moment regime occurs at elevated temperatures. Note that in 
the quantum critical region, the renormalization-group treat- 
ments give the same exponents of the SCR.^^' ^*^' The simi- 
lar idea of the SCR has been used in the microscopic the- 
ory referred as the two-particle self-consistent (TPSC),^"*'^^' 
where the charge and the spin fluctuations in the longitudinal 
channel are determined by enforcing the local sum rules to 
satisfy the Pauli principle. The double occupancy in the the- 
ory is expressed empirically in terms of the renormalized irre- 
ducible vertex used to describe the spin fluctuation, which in- 
troduces the self-consistent relation of the fluctuations. Once 
the two-particle fluctuations are obtained, the single-particle 
self-energy is computed by the way similar to the paramagnon 
theories. 

In order to go beyond the FLEX, Bickers and the col- 
laborators have argued the full consistency among one and 
two-particle levels using the parquet formalism. '"'•"'^"■'^' Note 
that the parquet formalism with the unrenormalized bare ir- 
reducible vertices is equivalent to the FLEX up to the two- 
particle level. ^*' They have approximated the irreducible ver- 
tices as the contact-type pseudo potentials at the representa- 
tive point in the momentum-frequency space in order to solve 
the parquet equation practically. The quantitative agreements 
have been demonstrated by the comparison with the quantum 
monte carlo (QMC) data for the intermediate repulsions. Nev- 
ertheless, there exists a methodological difficulty in choosing 
the representative point in the 4- vector space for general mod- 
els, which is usually unknown a priori. For the impurity An- 
derson model Janis and the co-worker have found the semi- 
analytic solution for the simplified parquet equation in the 
critical regime of the spin fluctuations.^^'''"' They have repro- 
duced the correct Kondo scale in the strong-coupling limit. 

These findings strongly encourage to construct a reasonable 
microscopic theory for strongly correlated electron systems 
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using the parquet formalism with appropriate modifications to 
the previous attempts. The parquet formaUsm is also suitable 
for the purpose that the two-particle self-consistency should 
be enforced prior to the single-particle one in view of the 
loop-expansion argument of the renormalization-group the- 
ory. Furthermore, the parquet equation is to be thermodynam- 
ically consistent, since it can be derived by functional deriva- 
tives of the generating free-energy functional with respect 
to the two-particle vertex functions."*'* It is not a conserv- 
ing approximation in the sense of Baym and Kadanoff.'^'^^' 
Needless to say, conserving approximations do not necessar- 
ily mean quantitatively reliable approximations. Indeed, the 
FLEX shows the considerable deviations from the finite-size 
QMC results even in the weak-coupling region. More- 
over, one should notice that any expression for the self-energy 
following the Baym-Kadanoff scheme is inconsistent with 
the exact self-energy formula of the Schwinger-Dyson equa- 
tion. '^•^"''"''' The both relations could be satisfied simultane- 
ously, only when the exact expression for the irreducible ver- 
tices would be found. "^^^ 

In this paper, we propose yet another self-consistent fluc- 
tuation (SCF) theory based on the parquet formalism, which 
will be explained in the next section. The formulation of the 
SCF theory will be given in §3, where the fully antisym- 
metric full vertex is constructed in terms of the two-particle 
fluctuations in all channels to ensure the Pauli principle. Re- 
garding the Pauli principle, the present theory is contrast to 
the TPSC, which satisfies it implicitly through the local sum 
rules. These fluctuations are determined self-consistently, and 
they are used to improve the one-particle self-energy in the 
same spirit of the TPSC. Note that the numerical cost of the 
present theory is much less than that of the FLEX. In §4, we 
demonstrate the applications to the impurity Anderson model 
and the Hubbard model on a square lattice as a benchmark 
of the present theory. In the last section, the relation to the 
SCR theory in the QCP region is discussed."*^* The present 
theory gives the same critical exponents and it can describe a 
crossover to non-universal region away from the QCP without 
introducing any phenomenological cutoffs as in the SCR. We 
also argue the possible improvements of the present theory, 
and we summarize the paper. In Appendix we summarize the 
derivation of the parquet equation for the single-band model 
in the case of a nonlocal general two-body interaction con- 
serving the z-component of the spins. 

2. Parquet Formalism 

In this section, we introduce briefly the parquet formalism 
and set up the notational convention through it. A detailed 
pedagogical review for the parquet formalism can be found in 
the literature.-'^' 

2. 1 Full vertex and crossing symmetry 

Let us first define an effective interaction with full vertex in 
the particle-particle (PP) and the particle-hole (PH) channels 
in the form, 

PP channel = ^r''(12; 34)c"'"(l)c''(2)c(4)c(3), 

PH channel = ^r(12; 34)c"'"(l)c(2)c"'"(4)c(3), (1) 



(a) 2: q-k 



4: q-k' (b) 2: k-q 4:k'-q 



1: k 



3:k' 



l:k 



3:k' 



Fig. 1 . The full vertices in (a) the PP channel and (b) the PH channel. The 
four-vector momenta are assigned to each legs as shown in the figure. The 
PP and the PH pairs carry the momentum q from right to left. 



where the number in the parenthesis represents a set of in- 
dices of the band (orbital), the spin, the spatial coordinate 
and the imaginary time. The superscript "p" indicates a quan- 
tity concerning the PP channel. The repeated indices that do 
not appear in the left-hand side are assumed to be summed 
hereafter We define the fourier and inverse fourier transfor- 
mations between the spatial coordinates and the four-vector 
momenta k - {k, ia>„), k' - {k', ia>'„), q - (q, ie,,,) as 

A(12; 34) = 2A(/t,/t';^)e'['^ '-('^-'?)-2-'^'-3+(*'-«)-4], 

kk'q 

A{k,k';q) = ^ A(12; 34)e-'[* i-(*-?)-2-*'.3+(«:'-?).4]^ ^2) 

234 

where w„, wj, are the fermionic Matsubara frequencies, while 
e„, is the bosonic one. We have assigned the four momenta at 
each legs of the vertex as A{k, k'; q) = A{k, k - q; k', k' - q) 
in the PH channel and AP(k, k';q) = AP{k, q - k;k' ,q - k') m 
the PP channel. The assignment of the variables for the full 
vertices are schematically shown in Fig. L 

The antisymmetric nature of the fermion fields gives the 
following relations for the full vertices, 

TP{D) = -TP{T), T{D) = -r(C), 

YP{D) = -Y{P) = Y{X), (3) 

where we have introduced sets of coordinates/momenta as 

D = (12; 34), (k,k';q), 

T = (12; 43), (k,q-k';q), 

C = (13; 24), (k,k-q;k-k'), 

P = (14; 32), (k,k';k + k' -q), 

X = il3- 42), (k,q-k';k-k'), (4) 

for notational convenience. These relations are called the 
crossing symmetry, which is the direct consequence of the 
Pauli principle. Note that frequently used theories such as the 
HF-RPA and FLEX violate the crossing symmetry. 

2.2 Bethe-Salpeter equation 

The full vertex can be decomposed into multiple scatter- 
ings of the PP or the PH pair, respectively, in terms of the 
irreducible vertices, Y^p{D) and r"(D). Namely, so-called the 
Bethe-Salpeter (BS) equation in the PH channel is given by 

r(D) = r°(D) + <D(D), <i)(D) = - [rVor] (£)), (5) 

where ;^'o(12; 34) = -G(13)G(42) is the lowest-order two- 
particle green's function. r*'(12;34) - r°(43;21), is irre- 
ducible with respect to a vertical cut of two green's function 
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Fig. 2. The Bethe-Salpeter equations in (a) the vertical PH channel, and (b) 
the PP channel. F" and F"'' are irreducible with respect to a vertical cut of 
two green's function lines. 



lines. Here the "matrix product" is introduced as 

[AB](D) = A(12;56)B(56;34), ^^Mk, p;q)B(p,k';q). 

p 

(6) 

The identity matrix is given by 1(D) = 5(13)^(24), da'- 
Similarly, the BS equation in the PP channel is given by 

= r°''(D) + *P(Z)), = [r°Vor^](£>), (7) 

where i/ro(12;34) = G(13)G(24). The factor 1/2 is to avoid 
double counting of twisted diagrams. Exchanging the indices 
1 and 2 in the BS equation, we can show that F"'^ also satisfies 
the crossing symmetry (3), i.e., r°P(D) = -r°P(T). Note that 
the irreducible vertex in one channel is reducible in another 
in general. The Feynman diagrams for the BS equations are 
given in Fig. 2. 

Using the full vertices, we express two-particle green's 
functions in the PP and the PH channels as 

= (r,c(l)c(2)c"'(4)c"'(3)) 

= MD) - \ [«Aor°V] (D) = MD) - \ [<AorVo] (£>), 

XiD) = (r,c(l)c'(2)c(4)c'(3)) 

= XoiD) - [xo^x] (D) = Xo(D) - hoTxo] (D). (8) 

The ordinary two-particle correlation functions are obtained 
by setting 1=2 and 3=4, i.e., xil) = Y^kk' x(k,k';q) and 
^iq) = Zm' Hk,k';q), where Y^k = ^llk 2« with A^o being 
the number of lattice sites. 

2.3 Parquet equation 

Although the equations given in the previous subsections 
are formally exact, it is necessary to introduce physically 
motivated approximations to make them practically solvable. 
However, a simple approximation to the irreducible vertices 
often results in full vertices violating the crossing symmetry. 
In order to obtain a solution that maintains the crossing sym- 
metry, it is adequate to use a set of exact relations between the 
irreducible vertices. This is the Parquet equation,-*^'^^* which 
is give by 

F°(D) = y(D) - 0(0 - ^'(P), 



r^PiD) = y''(D) + (D(X) - (B(f ). 



(9) 
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Fig. 3. The Parquet equations for the iireducible vertices in (a) the PH 
channel, and (b) the PP channel. By the BS equation with these irreducible 
vertices, the resultant full vertices satisfy the crossing symmetry. 



Here we have introduced the fully irreducible vertex, yiD), in 
all channels, which satisfies the crossing symmetry like (3). 
A choice of yiD) is ambiguous. In the simplest case, an anti- 
symmetrized bare interaction U (D) is used for y{D), and it is 
often called as "basic" parquet equation. 

Using the BS equations, we reexpress the full vertices as 

F(D) = y(D) + <D(D) - (D(C) - ^^(P), 

rP(D)^yPiD) + 'i'(D) + <:)iX)-0(P), (10) 

where the crossing symmetry is apparent. It should be em- 
phasized that resultant full vertices always satisfy the crossing 
symmetry, irrespective of the explicit expressions of O and 
The Feynman diagrams of the parquet equations are shown in 
Fig. 3. 

The BS and the parquet equations constitute a set of the 
self-consistent integral equations among the full and the ir- 
reducible vertices for an arbitrary given one-particle green's 
function G{k). 

2.4 Self energy in terms of full vertex 

Once we obtain the full vertex, the self-energy is calculated 
through the Schwinger-Dyson equation as shown schemati- 
cally in Fig. 4, 

2(12) = Shf(12) + 2(12), 

Ihf(12) = U''(13;24)G(43), 



1(12) = -i [UP^f/oFP] (13;24)G(43). 



(11) 



This equation together with the BS and the parquet equations 
provide a complete set of the self-consistent equations for cor- 
related fermion systems. 

2.5 Case of Hubbard-type interaction with SU(2) symmetry 
In the presence of SU(2) symmetry in spin space, the two- 
particle quantities can be classified by the total spin of the PH 
and the PP pairs, and the one-particle quantities are indepen- 
dent of spin indices (See, Appendix in detail). Namely, the 
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Fig. 4. The Schwinger-Dyson equation for tlie self-energy in terms of the 
full vertex. 



charge (even-parity) and the spin (odd-parity) vertices belong 
to 5 = and 5 = 1 PH (PP) channels, respectively. For later 
convenience, we quote the relevant expressions in the case 
of the Hubbard-type bare interaction, v(12) - U6(12) with 
SU(2) symmetry. 

BS equation 

T^(D)^Yl(D) + Q>^(D), (^^c,s). 

Of (D) = - 2 r»(^, p; q)xo(p; qW^ip, k'- q), 

p 

rf(£)) = r°(D) + »F^(Z)), {^^e,o), 

= -^rO(A:, p- q)Mp\ qW((P, k'; q), (12) 
where ;^o(^; q) - -G(k)G{k - q) and ij/Qik; q) - G(k)G(q - k). 
Two-particle correlation function 

Xdl) = Xoiq) - ^Xaik; q)T^^(k, k'; q)x^(k';q) 

kk' 

= Xoiq) - Y^Xoik; qW^ik, k'; q)xo(k'; q), 

kk' 

Mq) = Mq) - ^ X Mk; ?)r°(^, k'; q)ip^(k'- q) 
^ kk' 

= Mq) - ^ y Mk; qW((k, k'; q)Mk'; q), (13) 

^ kk' 

where ;i^o(?) = I,kXo(k; q) and iffo(q) = Yjk Mk; q)- Note that 
the (pair) density operators with 1 / V2 are used in the defini- 
tion of the correlation functions, i.e.. 



1 



1 



1 V 1 V 

Pe ^ -^C^{ia-)al3Cp, Po^ -^C^{lCr^(T)apCp. (14) 

Basic parquet equation (y = U) 

r°(£') = r - ^ [<E>. + 30,] (C) + \vve- 3*P„] {P\ 

r?(D) = -y - i [(D, - (D,] (C) -\['Ve+ {P), 

r»(D) = 2y + i [O, - 30,] iX) + \{^c- 30.,] (P), 



1 



1 



KiD) ^^[^c + O.J {X) - 2 [I*. + I*,] {P). (15) 

Here y is used instead of U in order to keep track of the fully 
irreducible vertex. 



Crossing symmetry 

Te{D) = r,(r), Yo(D) = -r„(r), 

TciD) = [rc + 3r,] (C), r,(D) = - ^ [r. - r,] (C), 

Tm = i [r, - 3r,] (P) = ^ [Fc - 3r,] (x), 

ToiD) = [r, + r,] (P) = ^ [r. + r,] (X). (i6) 



Self energy 



Zik) ^Un-^Y, 'f'o(p-^ + 3r„] (p, k; q)G(q-k), (17) 

where n = Yjk G(k)e'""'^+ is the electron density per spin. 

With this preliminary, we will discuss a practical approxi- 
mation of the parquet formalism in the next section. 

3. Self-Consistent Fluctuation Approximation 

The parquet formalism is a reliable nontrivial renormal- 
ization scheme, where two-particle self-consistency and the 
crossing symmetry are respected. However, a practical solu- 
tion is limited to very small systems since the central non- 
linear integral equations with three 4-vector variables are too 
large to handle even at modern numerical computation facili- 
ties.'*'** Therefore, a practical simplification preserving essen- 
tial ingredient of the parquet formalism is necessary. 

In order to simplify the parquet formalism, we introduce 
the contact-type effective irreducible vertices in the charge, 
the spin, the even-parity and the odd-parity channels, 

r;'(D)~A,, r?(Z))~-A„ 



rO(£)) ~ 2A„ rO(D) ~ 2A„ = 0, 



(18) 



which act as the renormalized vertices on an average, and 
they will be determined later without losing two-particle self- 
consistency and the crossing symmetry. The minus sign and 
the factor 2 are only for the notational simplicity of the fol- 
lowing equations. Note that the effective vertex in the odd- 
parity channel identically vanishes due to the crossing sym- 
metry, r"(Z)) - -r°(r). with this approximation, the formal 
solution of the BS equation (12) takes the form, 

<&f (D) ~ -Alxdq) = = c> 



(19) 



which are schematically shown in Fig. 5. Here, we have used 
the first expressions of (13) to yield the two-particle correla- 
tion functions, 

, , Xoiq) , . Xoiq) 
Xciq) = -r—K — TT' Xsiq) = 



^eiq) = 



1 + AcXoiq) 

^oiq) 



1 -A^xoiqY 



(20) 



1 + A,t//oiq) 

Using 'S>^iq) and ^eiq) and the BS equation (12), the full ver- 
tices from the parquet equation (15) read, 

T.iD) = y + <D,(^) - ^ [1>c + 30,] (0 + »P,(Q), 
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Fig. 5. The higher-order collective interactions in terms of the two-particle 
correlation functions. 



Fig. 6. The relation between the elfective iiTeducible vertices and the ver- 
tex coiTections. The (/-average of this relation gives the self-consistent 
equation to determine A. 



r,(D) = 2[r + + ^l^c- 3(S>,] (0 + ^ [O, - 3(DJ (fi), 



ToiD) = ^ [O. + O,] (Q) - i [O, + <D,] (Q), 



(21) 



where we have introduced two "bosonic" variables, Q - k-k' 
and Q - k + k' - q. The Pauli principle can be confirmed 
directly by these expressions. 

Substituting these vertices into the second definition of the 
two-particle green's functions (13), we obtain 

xM) = xoiq) -xaiqfly + ^M)} - ^xdq), 
Xsiq) = xo(q) +XQ(qf[7 - '^siq)] + Axs(q), 
Mq) = Mq) - Mqfiy + "Veiq)] - Ai/r,(^), 

^oiq) = ^oiq), (22) 
where the vertex corrections are defined as 

AXc(q) = ^ [AlMq) + 3AlMq)] - A^^g(q), 
^Xsiq) = -5 [Alfciq) - A]fM)\ - Alg(q), 

^^e(q) \AlhM) - 3A%(q)] . (23) 

Here we have defined the integrals corresponding to the Maki- 
Thompson processes as 



Mq) = Y^xoik; q)xd^ - k')xo(k';q), 

kk' 

g(q) = ^A'o(fc; q)i/^e(k + k' - q)XQ{k'\ q), 

kk' 

^dq) = X ^^^^^ ~ ^')'Ao(^'; q)- 



(24) 



kk' 



Equating (20) and (22), we obtain 

[Af - y\xQ{qf = A^f(^), 
[A, - y] <p^^{qf = Mie(q), 



(25) 



which are schematically shown in Fig. 6. Since we attempt to 
find a solution within the contact-type vertex approximation, 
the above equations make sense with an appropriate averaging 
procedure with respect to q. We adopt a simple average over 
q, and we finally obtain the self-consistent condition for the 
effective vertices, 

Y^qMleiq) 



Af = y H- — ^, K 



y + 



(26) 



Y^qXoiq) IjqMq) 

Note that the q average of f^(q), g(q) and h^(q) appeared in 




+ 




1> 



+ *- 



+ *- 



Fig. 7. The self-energy in the self-consistent fluctuation approximation. 
Each graphs represent the contributions from the 2nd-order, the longitu- 
dinal, the transverse and the PP fluctuations, respectively. 



the vertex corrections (23) reduces to the simple average, 

YjMi^ = Z'^f^^^ " T^^oiqfxdq)^ 



Y^g(q)^Y^iPo{qf<PM)- 



(27) 



By solving (23), (26) and (27) self-consistently for given 
G(k), we obtain the renormalized vertices, A^, A, and A^,, and 
simultaneously the two-particle correlation functions and the 
fully antisymmetric vertices via (21). 

Substituting the full vertices (21) into (17) and transform- 
ing the variables, we obtain the self-energy 

nk) = un + u Ylxo(q) {r + ^ [1>c - 3<D,] (^)| G(k - q) 



- Mq)'^e{q)G(q - k)\. (28) 

The self-energy in this approximation is schematically shown 
in Fig. 7. Note that when we neglect the vertex corrections, 
A;^fc = ^Xs - M'e - 0, the present formalism reduces to the 
FLEX. 

Let us now comment on the self-consistent procedure be- 
tween two-particle and one-particle levels. In the two-particle 
self-consistency, the one-particle propagator is an external in- 
put. It is well known that the one-particle self-consistency 
without using the consistent frequency-dependent irreducible 
vertices always results in a universal high-frequency behav- 
ior and the Hubbard satellite bands are completely smeared 
out.^"^'^^' This is because the high-frequency asymptotic be- 
havior sets in, not above the bandwidth but above an incorrect 
energy scale, U, which is extensively discussed by Vilk and 
Tremblay,-'"^' and Janis and Augustinsky.-*^^ 

In order to balance the use of the contact-type irreducible 
vertices, A^ and A^,, we use the HE propagator as an input 
for the two-particle self-consistency, and the one-particle self- 
energy is updated only once. Specifically, we first determine 
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the HF propagator Goik) by adjusting the chemical potential 
to match the given electron density per spin, n. Using Goik), 
we determine , As and A^, by solving the two-particle self- 
consistent equations. Then, the self-energy l,(k) is updated via 
the Schwinger-Dyson equation, and the renormalized chemi- 
cal potential, d/j, is determined-''*^ by Yik G(k)e"^"^* - n with 
G(fc)-i = Go(/t)-i - Uk) + 6fi. 

Empirically, the contribution from the PP fluctuation seems 
to be slightly overestimated, and it yields unphysical behav- 
iors in one-particle spectrum. By neglecting the last term in 
(28), the SCF works quite well as will be shown in the next 
section. 

4. Application to Fundamental Models 

In this section, we consider two fundamental models of 
many-body fermion systems, i.e., the impurity Anderson 
model and the Hubbard model on a square lattice, as applica- 
tions of the SCF approximation. In the latter, the first Brillouin 
zone is discretized by 128 x 128 meshes, and 2048 Matsub- 
ara frequencies are used in both models. To obtain the real- 
frequency spectrum, we adopt the Pade interpolation to carry 
out an analytic continuation. 

4.1 Impurity Anderson model 

The impurity Anderson model is given by 



ka- 



+ ^ Z i'lj- + flc,^) + Unnrifi, (29) 



where the / level is measured from -f//2, and n/a- - fa-fa-- 
We use the practically infinite bandwidth of the conduction 
electron, so that the Hartree green's function is given by 



Go(/w„) 



f/(nH- l/2)-H/Ao, 



where n^ - T Ym G{ia)„)e"^"''^+ . The hybridization strength is 
Ao = npoV^, which is taken as a unit of energy. 

Let us consider the symmetric case {Ef - 0), in which 
GoiicJn) and G{ia>„) are pure imaginary, and « = «h = 1/2, 
Xoiiem) = Mi^m) and A,. = A^ hold. 

The U dependence of the eff'ective irreducible vertices at 
T - O.OlAo are shown in Fig. 8. The charge (even-parity) 
vertex is roughly the order of U, while the spin vertex is 
strongly renormalized as U increases. In the HF-RPA theory, 
the Stoner condition is given by U/ttAq = 1, while the spin 
vertex in the present theory never reaches the magnetic insta- 
bility condition 6=1- AsXoi^) = as it should in the essen- 
tially zero-dimensional system. Note that 6 gives a measure 
of the Kondo temperature, i.e., 6 ~ xo(0)Tk. 

The T dependences of A^. and A, at f/ = AnAo are shown in 
Fig. 9. Both the charge and the spin vertices have very weak 
temperature dependence, and the T dependence of the suscep- 
tibilities thus mainly comes from that of ;^o('fm)- 

The T dependence of the static charge and spin suscepti- 
bilities, ;^^c(0) and;t's(0) at t/ = 47rAo are shown in Fig. 10. 
The charge fluctuation is suppressed below T ~ U, while the 
saturation of the spin susceptibility occurs below T ~ 0.1 Ao 
reflecting the singlet formation of the Kondo effect. The exact 
expression of the Kondo temperature in Kondo regime 
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Fig. 8. (Color online) The U dependence of A^. and A, for T/Aq = 0.01. 
Note tiiat the PP vertex A^ is identical with A, for the symmetric case. 
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Fig. 9. (Color online) The T dependence of A^ and A^ at U/kAq = 4. 



(U, \Ef\ — > oo) is given by 



Tk 
Ao 



71 U 

ItiAq 



exp 



' 8 ;rAo 



1 - 



2£/ 
U 



For the intermediate coupling, U - AttAq, Tk/Ao ^ 0.018 and 
the present theory slightly overestimates the Kondo tempera- 
ture. 

The dynamical spin susceptibility is shown in Fig. 11. As 
the temperature decreases, the intensity of Xs(£) transfers to 
the lower-energy region. The low-energy limit of Im;^'j(e) is 
proportional to e as the FL theory. 

Let us now turn to the one-particle density of states, p{a)) - 
-lmG{cij)/n. The U dependence of p{aj) at T = 0.0 lAo is 
shown in Fig. 12. The central Kondo quasiparticle peak is well 
developed and its width becomes narrower as U increases. 
The position of the Hubbard satellite peaks is slightly smaller 
than the values of the atomic limit, +U/2. The overall features 
of the DOS are well reproduced and is quantitatively similar 
to the best know results by the numerical renormalization- 
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Fig. 10. (Color online) The r dependence ot>^(0)and^,(0) at f//Mo = 4. Fig l2. (Color online) The J/ dependence of p(w) = -\mG(a>)ln al 

T/Ao = 0.01. 
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Fig. 11. (Color online) The r dependence of Im^,(e) at {//;rAo = 4. 

group method. 

The T dependence of p(co) atU - AjtAq is shown in Fig. 13. 
The Kondo peak begins to develop below Tk where the static 
spin susceptibility also begins to saturate. Note that the zero- 
temperature limit of p(0) is fixed as p(0) = 1 /ttAq due to the 
Friedel sum rule.^^^ 

The quantitative success of the density of states has a com- 
mon origin of the second-order perturbation theory in the 
symmetric case.*^"^^^ Namely, the dominant contribution of 
the fluctuations in (28) is the second-order term (the first term 
in the curly bracket). By the same reason, the overall features 
of p(a») in the asymmetric case (not shown here) are similar 
to those of the second-order perturbation theory.^'-^^' Nev- 
ertheless, the low-order perturbation theory is inapplicable 
in higher-dimensional systems in principle, where the phase 
transitions are possible to occur The present approach can 
describe potential phase transitions in higher dimensions, and 
at the same time the artificial instabilities of the mean-field 
theory can be successfully eliminated by the self-consistent 
treatment of the fluctuations. The concrete example of the 



Fig. 13. (Color online) The T dependence of p{oj) = -lmG{oj)/n for 
U/nAo = 4. 

two-dimensional Hubbard model will be shown in the next 
subsection. 

4.2 Hubbard model on square lattice 

Next, we consider the Hubbard model with the nearest- 
neighbor hopping, 

= ^ ekcl^Cko- + nijriii, (30) 

kcr i 

where eu - -2f[cos(^ta) -i- cosik^a)}, with a being the lattice 
constant. The number operator at / site with the spin cr is given 
by n,o- = '^I'o-'-Hr- following we use a = f = 1. In the 

particle-hole symmetric case at the half filling, we have the 
relations, G{k) = -G{k + Q), I.(k) = -I(fc + g), xuiq) = 
iAo(^ + Q) and Ac = Ap, where Q = (Q; 0) and Q - (n, n). 

Let us discuss the symmetric case at the half filling, n - 
1/2. Figures 14 and 15 show the U and T dependences of 
Ac and As- The overall behaviors are similar to those of the 
impurity Anderson model, e.g., Ac is of the order of U, As is 
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Fig. 14. (Color online) The U dependence of the elfective irreducible ver- 
tices for T = 0.2f. 
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Fig. 15. (Color online) The T dependence of the effective in'educible ver- 
tices for U = 8f. 



renormalized to eliminate the magnetic instability, and both 
of Ac and A, have weak T dependences. 

For T < 0.2t in Fig. 15, AsXoiQ,0) approaches to unity, 
but it never exceeds the value of the magnetic instability to- 
ward r = 0. In order to elucidate the Mermin-Wagner theo- 
rem, the T dependence of the inverse staggered susceptibility 
is shown in Fig. 16. l/Xs(Q,0) follows the Curie- Weiss be- 
havior in higher temperatures, and it begins to deviate below 
T = 0.2f. For later convenience, we define the fictitious Neel 
temperature, by extrapolating the Cruie-Weiss T depen- 
dence, corresponds to "mean-field" Neel temperature, and 
it is regarded as the true critical temperature in the presence of 
a weak three dimensionality. The inverse susceptibility is then 
approaching to zero with the exponential T dependence. The 
dashed lines in the figure show the fitting by T exp(-2b/T) 
or exp(-2b' /T) with or without the quantum correction.-" 
The relation between the present theory and the SCR in the 
quantum critical region will be discussed in the later section. 

The intensity map of Imxsici, e) for U - St and T = 0.18? 




0.25 



Fig. 16. (Color onhne) The T dependence of the inverse static AF spin sus- 
ceptibility, 1/^1(6,0). The dashed lines are the predicted T dependences 
■pg-lb/T g~2(> IT ^jjjj yyjjjjout the quantum correction. The inset is the 
logarithmic behavior of 1 IxAQ, 0) as a function of l/T. 
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Fig. 17. The intensity map for ImXsiq, f) along the high-symmetry line 
of the Brillouin zone, where X = {tt, 0), r(0, 0), M = (tt, tt) and M' = 
(7T/2,n/2). The critically enhanced paramagnon excitation is prominent 
with the sharp peak at M point. 



along the high-symmetry line of the Brillouin zone is shown 
in Fig. 17, where X = (7r,0), F = (0,0), M = (n,7T) and 
M' = {71 12, Ti 12). Since the spin susceptibility is of the RPA- 
type with the renormalized vertex A 5, the support of the 
continuum excitations is given by the possible combinations 
of the particle-hole pair excitations. Note that the critically 
enhanced collective paramagnon excitation is prominent at 
lower temperatures. 

The imaginary part of the dynamical staggered suscepti- 
bility for U - 8f is shown in Fig 18. It shows the strong T 
dependence to enhance the low-energy excitations as T de- 
creases. As will be shown below, the critical AF spin fluctua- 
tions result in the pseudogap behavior in the one-particle ex- 
citation, which should give a feedback to the low-energy exci- 
tation of hnxiQ, e) as well. Since the present theory lacks the 
self-consistency between the one-particle and the two-particle 
levels, no suppressions in the low-energy excitations appears 
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Fig. 19. (Color online) The T dependence of np{cS) = - 2* '^)- 

as a feedback. The full self-consistent treatment remains an 
open and very challenging problem. 

Let us now turn to the one-particle properties. The T de- 
pendence of the DOS for U = 8? is shown in Fig. 19. As T 
decreases the central quasiparticle peak develops, and simul- 
taneously the pseudogap in the quasiparticle peak becomes 
deeper at w = 0. Note that the development of the pseudo- 
gap begins at higher temperature where the staggered spin 
susceptibility deviates from the Curie-Weiss behavior. The 
critical AF spin fluctuations have little influence on the inco- 
herent part, and the transfer of the weight mainly takes place 
within the quasiparticle peak. This tendency is in contrast to 
that of the DMFT with nonlocal corrections where consider- 
able amount of the weight transfers to the incoherent Hubbard 
satellite peaks from the quasiparticle peak.'^-^^'^^' 

The spectral intensity map of A(k, uS) - -Im G(k, uS) along 
the high-symmetry line of the Brillouin zone for f/ = 8f at 
r = O.lSf is shown in Fig. 20. Since there is a perfect nesting 
at half filling and the magnetic Brillouin zone coincides with 
the Fermi surface, the pseudogap behaviors appear in the co- 
herent quasiparticle bands everywhere on the Fermi surface 
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Fig. 20. (Color online) The A{k,(ii) = -\mG(k,tij) along the high- 
symmetry line of the Brillouin zone. 

(see, the X-M' line). The overall structure of A(k, a)) is very 
similar to that obtained by the DMFT with nonlocal correc- 
tions.^^- 

Now, let us consider the hole doping and the associated su- 
perconducting transitions. In the present theory, the even-pair 
and the odd-pair irreducible vertices are given by 

r«(D) = 2r + i [0,(0 - 30,(0] + i [cD,(Q) - 30,(Q)] , 
r°(D) = i [0,(e) + 0,(2)] - i [0,(Q) + (5,(Q)] . 

(31) 

Provided the irreducible vertices, we determine the supercon- 
ducting transition temperature Tc by the linearized BS equa- 
tion, 

Acp^ik) --\Yj ^)\Gik'fMk'\ = e, o), (32) 

p 

where A - \ signals the superconducting instability of the 
normal state, and its eigenvector (p({k) provides the symme- 
try of the gap function.'^' Here, Oc(g') and ^siq) are given by 
(19). For the approximate vertices (32), we easily solve the 
linearized BS equation by using the power method with the 
fast fourier transformation (FFT) algorithm. On the Hubbard 
model, we only consider the t/^:_,,2-type (Biij-type) order pa- 
rameter, which is the favorable symmetry mediated by the AF 
spin fluctuations. 

Away from the half filling, the "mean-field" rj*| decreases 
and it vanishes at n ~ 0.43 as shown in Fig. 22. Solving 
the linearized BS equation, the maximum eigenvalue A does 
not reach unity at the lowest temperature (not shown) for the 
doping range, 0.3 < n < 0.5. As compared with the FLEX, 
which is often used to discuss the superconducting instabil- 
ity, the present theory uses the much smaller coupling con- 
stant to the spin fluctuations. A, renormalized from U. It is 
the main cause to suppress Tc in the present work. However, 
the constant approximation for the irreducible vertices may 
underestimate the coupling constant to the spin fluctuations. 
If we considered the momentum and/or the frequency depen- 
dence, Asiq) would have a minimum at (Q; 0) and it would 
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Fig. 21. (Color online) The T dependence of the eigenvalue of the lin- 
earized BS equation. The dj.2_j,2-type of the order parameter is assumed. 
In calculating Tc, the spin vertex is multiplied by a factor g = 2.5, i.e., 
A J gKs (see text in detail). 



approach to U away from (Q; 0). Since the strongest fluctua- 
tion tends to break the Cooper pairs, the electrons make use 
of a broad range of fluctuations around the maximum point of 
Xsiq) where the coupUng constant hsiq) would be larger than 
its averaged value. In order to mimic this effect, we merely 
replace A, with gK^ {g > 1). The simple enhancement of A, 
of course needs to be justified by more elaborate calculations 
with ^-dependent vertices. Nevertheless, it is useful to grasp 
the doping dependence of Tc- For instance, setting g = 2.5 
(gAs ~ Q.6U), we have the T dependence of the eigenvalue 
as shown in Fig. 21. In this case, the eigenvalues reach unity 
with the non-monotonous doping dependence. Note that Tc is 
suppressed strongly when the system goes into the pseudo- 
gap region. The superconducting phase boundary is indicated 
by the closed triangle in Fig. 22. In order to elucidate the re- 
tardation effect in the self-energy, we also calculate Tc using 
the bare propagator, i.e., G Go, in the linearized BS ker- 
nel (32) with no enhancement factor, g = 1 . In contrast to the 
decreasing behavior of Tc to the pseudogap region, we obtain 
the monotonously increasing Tc as n increases (the open tri- 
angle in Fig. 22). Therefore, we conclude that the suppression 
of Tc is mainly due to the retardation effect of the self-energy. 

5. Discussions and Summary 

5.1 Relation to the SCR theory in the vicinity of the QCP 

We discuss the relation between the SCR and the present 
theory near the QCP. The core of the SCR theory is the self- 
consistent equation, 

5{T)^5o + g(x.M))T, (33) 
where the dynamical spin susceptibility is approximated as 



Xs{Q* + qJe„r)-' ~6iT)+A\q\^ + C- 



z-2 ' 



(34) 



in the critical region, with z being the dynamical exponent. 
The Q* is a magnetic ordering vector. The bracket means the 
g'-average, (■ ■ ■)7- = (T/Nq) 2m 2? ' ' s and the subscript T 
represents that the significant T dependence is included in the 
average. The coefficient g is related to the mode-mode cou- 
pling constant in the 4th-order free-energy functional. S(T) - 
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Fig. 22. (Color online) The T-n phase diagram. The is the crossover 
temperature below which the pseudogap develops. In calculating T^, the 
spin vertex is multiplied by a factor g, i.e.. A, ~> g^s, and using either G 
or Go (see text in detail). 



XsiQ*,Oy^ is a square of the inverse correlation length, and 
it is regarded as a measure of the distance from the magnetic 
QCR The peculiar T dependence of 6{T) determined by (33) 
governs the critical exponents of various quantities in the QCP 
region. Setting 5 = Q atT - 0, we obtain the condition for do 
to give the QCP, 



6q 



r (xs(q)}T=o ■ 



(35) 



In the present theory, if we have the single peak in x.^{ci, 0) 
at ^ = g*, we have the relation. 



6{T) = Xo(Q\Or' - As(T), 



(36) 



where the weak T dependence in xo can be neglected in the 
QCP region as long as additional singularities such as the van 
Hove singularity and the perfect nesting are not involved. The 
self-consistent equation (26) for A, in the critical region may 
read 



A, = f/ - c -H ffAj (xM)}t ■■ 



(37) 



where c is constant contributions from the charge and the PR 
fluctuations, and a is an (9(1) constant. In the leading order in 
(5, it reduces to 



6{T)^5'^ + g' (xM))t^ 



(38) 



which is the same form as (33) with different coefficients. The 
formal equivalence of the self-consistent equations thus leads 
to the same exponent of 6{T), and other relevant exponents in 
the QCP region. 

Note that in general the basic parquet equations lead to the 
same critical exponents of a self-consistent 1 jN expansion for 
the corresponding order-parameter field theories where is 
the number of order-parameter components.^^^ In contrast to a 
theory that handles a mode divergence in a single channel, the 
present theory naturally takes account of interference effects 
among different modes in the presence of the plural diver- 
gences in different channels or multiple characteristic wave 
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vector g*. 



and it is a very challenging open problem. 



5.2 Possible improvements for the present theory 

Let us consider possible improvements for the effective ir- 
reducible vertices. The extension of ^-dependent vertex, A(q) 
is formally straightforward. Namely, the self-consistent equa- 
tions (26) now become 

Af(«) =7+ ^z7Tr= --J, A,(q) =y+ 

(39) 

where the q average has been replaced by the €„ average. Al- 
though a formal change is very little, it increases a practi- 
cal computational effort considerably in the evaluation of the 
vertex correction terms (23) where the simple FFT algorithm 
cannot be apphed.^"^^ 

One may consider that the extension of the full ^'-dependent 
A(^) would also be straightforward in a similar way. However, 
the self-consistent equations. 



K^{q) = y + 



Kiq) = 7 + 



Ai//e(q) 



xoiq)^ ' ' foiqf ' 

have an apparent inconsistency. For example, in the high- 
frequency limit, the left-hand sides should reduce to the bare 
interaction, while the right-hand sides diverge as e^. This may 
be due to the fact that only the collective fluctuations do not 
play a dominant role in the high frequencies, but various in- 
dividual excitations equally play a role. All these fluctuations 
tend to cancel with each other yielding the recovery of the 
effective irreducible vertices to the bare interaction. There- 
fore, a construction of the frequency-dependent vertex theory 
in this direction seems to be a very challenging problem, un- 
less an appropriate parametrization for the irreducible vertices 
is found. 

5.3 Summary 

We have proposed the self-consistent fluctuation theory that 
treats the self-consistency and the antisymmetric nature of 
the two-particle fluctuations microscopically and appropri- 
ately. The applications to the impurity Anderson model and 
the Hubbard model on a square lattice have shown the self- 
consistent treatment of the fluctuations successfully elimi- 
nates the finite-temperature magnetic instability of theories 
without vertex corrections, and it gives the natural energy 
scale of the spin fluctuations. In the quantum critical region, 
the present theory is shown to give the same critical exponents 
as the SCR theory where a single mode diverges in the spin 
channel. 

The self-energy and the corresponding one-particle spec- 
trum reflect the developed fluctuations in a consistent way 
via the Schwinger-Dyson equation, and it can be described 
marginally the Hubbard satelhte peaks due to the local spin 
fluctuations. Since the DOS at the fermi energy and the dou- 
ble occupancy as a function of the temperature or of the inter- 
action strength do not show a discontinuous change, but show 
a smooth crossover, the Mott metal-insulator transition can- 
not be described in the present theory. Regarding this point, 
it cannot describe the energy scale of J ~ t^/U. The present 
theory is thus limited to use in the range of the intermedi- 
ate coupling, U < W.To discuss the strong coupling regime, 
an extension to the frequency-dependent vertices is necessary. 
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Appendix: Parquet Equation for Single-Band Model 
with Nonlocal Two-Body Interaction 

In this appendix, we consider the explicit spin dependences 
of the parquet formalism for single-band models. For this 
purpose, we explicitly write the spin indices as 1 — > (la), 
G(12) ^ Gap(l2), r(12;34) ^ rap;ys(12;34) and so on, 
where the greek letters represent the spin indices. Let us con- 
sider the two-body interaction that conserves the z-component 
of the spins, 



1 



Hm = -v„^(12)4(l)c;(2)c^(2)c„(l). (A-1) 

The anti- symmetrized interactions in the PH and the PP chan- 
nels are expressed as 

f/„^;^5(12;34) = v„5(14) [6(12^)6(3^4)6^^6^6 

-6(l3^)6(2^4)6ay6i}s] , 

t^«ftr5(12; 34) = v„^(12) [6(13M24^)6ay6ps 

-6(U^)S(23^)6as6i3y]iA-2) 

where the subscript of the coordinate, e.g. 1+, specifies the 
limit of the imaginary time as ri + 0. 

When a system conserves the z-component of the spins, 
the one-particle quantities are diagonal in the spin space, e.g.. 
Gap = GaSajs- For two-particlc quantities, we can decompose 
the spin dependences into the "direct" and "exchange" parts 
only within the spin space as 

Kfi:y6(D) - [K'^'s,ySp6+APj;X,6i,y](D), 

Aafi-ys(D) = [A'^^^6afi6ys + A^^6^6^s\ (D). (A-3) 

Here, A (A'') represents any of two-particle quantities such 
as the full vertex F (F''), the matrix product O (*P) and so 
on. Note that the bare two-particle green's functions are the 
special cases with A'^^^ = and A'-^ = 0. 

The explicit spin dependences of the parquet equation (9) 
are thus obtained by replacing 

(D(C) ^ OarMO = [^al^a06ys + <1>^^^<5„^<5^,1 (C), 

(D(X) ^ ^ay-M^) = i^ap^ay^e^ + ^^^^asSpy] (X), 
a>(P) ^ ^aS-yp(P) = [<l>L^^«rV + ^L*^ 

ad^l3y\(P)- (A4) 

The relevant matrix products are given by 
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3>S(^) = -KW.r«l(D), 

(A-5) 

Similarly, we obtain the self-energy by using (A-2) as, 

PI 

(A-6) 

where n„ - Yjk Ga{k)e'""^* is the number density for spin a. 

Now, let us consider the case of SU(2) symmetry in spin 
space. In this case, v^p, A^^'^^ and A^^^ do not depend on the 
spin indices. Using the identity, 

= ^ \Sal36sy + O-afi ■ 0"5y] , (A-7) 

and the transpose property, (icr^(T)i2 = (jcr^cr)2i, {i<T^)i2 = 
-(icr^)2i, (A-3) can be expressed as 

(A-8) 

where we have defined 

4 _ apW _ apM _ aP _ aP 

4 _ apW . apM — aP — aP + aP 

Ac = lA^^ + A« = Att;rT + Anai, 

A. = A^-'"' = At|;t^ = An;n - A„;^. (A-9) 

There is no ambiguity to omit the superscript "p" in the case 
of SU(2) symmetry. It is obvious from (A-8) that A,,, A^, Ac 
and A , represent the two-particle quantities in the even-parity, 
the odd-parity, the charge and the spin channels, respectively. 

Inverting these relations, and substituting them into (A-5), 
we obtain, 

m^iP) = [r^^oFfl {D\ (f = e, o). (A-10) 

Similarly, we have the parquet equations in the case of SU(2) 
syimnetry, 

T^(Z)) = jciB) - ^ [3>c + 3<D,] (C) + ^ ["P. - 3Y„] (P), 
r^(£)) = y,(U) -\\<^c- O,] (C) - ^ + (Z'), 
I^(D) = r,(D) + ^ [3>c - 30,] (X) + ^ [O, - 30,] (P), 



r^(D) = y„(D) + ^ [3>c + 3>,] (X) - ^ [O, + OJ (P). 

(A-11) 

In the basic parquet equations, i.e., jiD) = U (D), we obtain 

yc{D) = 2v{q)-vik-k'), 

yAD) = -vik-k'), 

jeiD) = vik -k') + v{k + k' - q), 

yo{D) = v{k-k!)-v{k + k! -q). (A-12) 

The expression of the self-energy in the case of SU(2) is 
given by 

2hf(A:) = v(0)n - [^(^) - 1^(0)] G{k - q), 

^(k) = 2 ''^^ " 'P^^o(p; q) [r. + 3rj (p, t, q)G{q - k). 
PI 

(A- 13) 
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